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Abstract 

> . 

, The square root of the positive definite hermitian operator DJ^Dyj in Neuberger's 

proposal of exactly massless quarks on the lattice is implemented by the recursion for- 
mula Y k+1 = \{Y k + DlD w Y^ 1 ) with y = I, where Y k 2 converges to Dj^Dyj quadrat- 
qq ■ ically. The spectrum of the lattice Dirac operator for single massless fermion in two 

dimensional background U(l) gauge fields is investigated. For smooth background 
+^ ■ gauge fields with non-zero topological charge, the exact zero modes with definite chi- 

'-h , rality are reproduced to a very high precision on a finite lattice and the Index Theorem 

Q_i! is satisfied exactly. The fermionic determinants are also computed and they are in good 

agreement with the continuum exact solution. 

> : 

^ ! 1 Introduction 

p.: 

In recent publications |2|], Neuberger suggested that in a vector-like gauge theory the 
lattice Dirac operator D for exactly massless quarks can be represented by a finite matrix of 
fixed shape, without undesired doubling and no need for any fine tuning. The lattice Dirac 
fermion operator of Neuberger's proposal of exactly massless fermion is 

D h = l + V, V = D w {DlD w )- 1 ' 2 (1) 

where D w denotes the standard Wilson-Dirac lattice fermion operator with negative mass 
term, and V is unitary ( = V~ x ). Dh was derived based on the observation that the 
overlap 0, [|, |5| for Dirac fermion in odd dimensions || can be written as the determinant 
of a finite matrix of fixed shape. The most remarkable feature of is that it satisfies the 
Ginsparg- Wilson relation 

D 75 + l5 D = D l5 D (2) 

and Dh is so far the only known explicit solution of eq.(||). Ginsparg- Wilson relation was 
derived in 1981 as the remnant of chiral symmetry on the lattice after blocking a chirally 



1 



symmetric theory with a chirality breaking local renormalization group transformation. The 
original Ginsparg- Wilson relation is in fact more general than eq.(|2|) and constitutes a matrix 
R which is local in the position space but diagonal in the Dirac space 

^75 + I5D = D l5 RD (3) 

It serves as the criterion for breaking the continuum chiral symmetry on the lattice while 
preserving the exact masslessness and the continuum axial anomaly. Recently, Hasenfratz, 
Laliena and Niedermayer [pj] , and Liischer Jl2| explicitly showed that any D satisfying 



the Ginsparg- Wilson relation plus some reasonable assumptions such as locality and free of 



species doubling must obey the Index Theorem on the lattice. Furthermore, Liischer |L2 
discovered that any D satisfying Ginsparg- Wilson relation implies an exact symmetry of 
the fermion action which may be regarded as a lattice form of chiral symmetry, reproduces 
the correct anomaly and the Index Theorem on the lattice. The Nielsen-Ninomiya theorem 
Tj| can be circumvented if the continuum chiral symmetry of the fermion is replaced by 



the Ginsparg- Wilson relation which is the chiral symmetry realized on the lattice. Recently, 
Narayanan [14] has shown that the fermionic determinant of Dirac fermion operator D in 



the form D = 1 + V and satisfying the Ginsparg- Wilson relation can be factorized into two 
factors which are complex conjugate of each other, corresponding to those of left-handed 
and right-handed Weyl fermions. This could imply that the dynamical gauge theory of 
single massless Dirac fermion using is amenable to Hybrid Monte Carlo simulations. 
However, the most challenging part of Neuberger's proposal is the implementation of the 
square root of the positive defmtite hermitian operator D^Du,. The approximation of Dh 
has been studied by Neuberger |T5[], but it sacrifices strictly masslessness. In this paper, we 
attempt to implement the square root operation on the positive definite hermition operator 
Dl u D w by the following recursion formula JTB] 



n+i = \{Y k + D^DJTj- 1 ), Y = 1 (4) 

It can be shown that Y" fc 2 converges to D^D^, quadratically. The main purpose of this paper 
is to investigate the spectrum of Dh to see to what extent this implementation can reproduce 
the exact zero modes with definite chirality as well as the realization of the Index Theorem 
on a finite lattice. The fermionic determinants are also computed and compared with the 
continuum exact solutions. 



2 Massless Fermion Action 

Neuberger's lattice fermion operator for exactly massless fermion has been given in eq. ([!]). 
The negative mass term in the standard Wilson-Dirac lattice fermion operator D w can be 
chosen to be any value in the range (—1, 0). A different value of m corresponds to a different 
renormalization for the observables. In the following, we shall restrict our discussions to the 
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case of m — — 1. Then the Wilson-Dirac operator becomes 

1 



£>„ = -! + -[7^(V* + V„) - V* VJ (5) 



where and are the forward and backward difference operators defined in the following, 

v^(x) = u^(x)ip(x + A) - i>( x ) 

V*^(x) = ip(x) - Ufa - p)tp(x - p) 
Then the lattice action of single massless fermion in background gauge field is 

A h = J2J2^(x)D h (x,y)ij(y) (6) 

x y 

where the Dirac indices are suppressed. 

On a torus ( x^ e [0,L M ],/i = 1,2 ), the U(l) gauge fields can be decomposed into 
global, harmonic and local parts. In this paper we use the following decomposition. 

AW = -^ + ^ + 4 0, sm(^x 2 ) (7) 

L\L 2 L\ L 2 

A2{x) = ^i + A f A< ^pi Xl) (8) 

L 2 L, 1 
where the global part is characterized by the topological charge 

Q = 7 kJ £xFl2 (9) 

which must be an integer. The harmonic parts are parameterized by two constants hi and 
h 2 . The local parts are chosen to be sinusoidal fluctuatations with amplitudes and A 2 °^ 
and frequencies and ^j^- where n\ and n 2 are integers. The discontinuity of A\(x) at 
X2 = L 2 due to the global part only amounts to a gauge transformation. The field strength 
F12 = d\A 2 — d 2 A 1 is continuous on the torus. To transcribe the continuum gauge fields to 
the lattice, we take the lattice sites at x M = 0, a, (N^ — l)a, where a is the lattice spacing 
and L M = N^a is the lattice size. Then the link variables are 

U x {x) = exp[L4i(x)a] (10) 

2nQxi 

U 2 (x) = exp[L4 2 (x)a + i — 4 2i(A r 2 _i )a ] (11) 

The last term in the exponent of U 2 (x) is included to ensure that the field strength F\ 2 which 
is defined by the ordered product of link variables around a plaquette is continuous on the 
torus. 

The fermion propagator Sf(x, y) is defined by 

Sf(x, y) = ^J]J dij(z)d^(z) e~ A » if>(x)$(y) (12) 
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where 

Z = f Y[diP(z)diJj(z)e- Ah (13) 



In background gauge fields of zero topological charge ((5 = 0), the fermion propagator is 

S F (x,y) = D h \x,y) (14) 

The free fermion propagator in momentum space is 

c~J°) (n) - ° i _ i0 E M 7„sin(;M 

{p) ~2 I ~ m 2[N(p) + u(p)} (15) 

where a is the lattice spacing and 

u(p) = 1-J3 cos (Pn a ) ( 16 ) 



JV(p) = ^w 2 (p)+Esin 2 (p, 1 a) (17) 

The constant term |1 which vanishes in the continuum limit is expected to appear in eq. (p~5|) 
such that D^ 1 satisfies the Ginsparg- Wilson relation 

75/T 1 + D~ x lf> = a 75 (18) 

which is equivalent to Eq.(||]). The denominator 2[N(p) + u(p)] in eq . (|15D has only one zero 
at pi = P2 = for the entire Brillouin zone and its expansion around pi = p^ = is 

2[N(p) + u{p)} = a 2 {p\ + pi) + O(oV) (19) 

Therefore the free fermion propagator is free of doublers and has the correct continuum limit. 
We remark that it is trivial to construct D^ 1 satisfying Ginsparg- Wilson relation. Any D^ 1 
in the form 

D- 1 (x,y) = ^l + J2l,S^y) (20) 

must satisfy eq. (|18"D . However, the additional requirements that the resulting lattice fermion 
action (D^ 1 )^ 1 is local, free of species doubling and satisfying the Index Theorem would make 
the task rather non-trivial. So far, is the only known explicit solution fulfilling all these 
requirements. 



3 Zero Modes and The Index Theorem 

In the continuum the Dirac operator of massless fermions in a smooth background gauge field 
with non-zero topological charge Q has zero eigenvalues and the corresponding eigenfunctions 
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are chiral. The Index Theorem || asserts that the difference of the number of left-handed and 
right-handed zero modes is equal to the topological charge of the gauge field configuration : 

n_-n+ = Q (21) 

In two dimensions, the so called Vanishing Theorem [[J also holds : 

Q > =>- n + = and Q = 

Q < n_ = and Q = -n+ (22) 



Recently, Hasenfratz, Laliena and Niedermayer |TTJ, and and Liischer ]T2[ explicitly showed 



that any lattice Dirac fermion operator D satisfying the Ginsparg- Wilson relation eq.(0) 
must obey the Index Theorem exactly. Since is an explicit solution of Ginsparg- Wilson 
relation, it must obey the Index Theorem exactly. However, on a finite lattice, after we 
implement the inverse square root operator in by the recursion formula eq.(^), it is not 
known a priori how the zero modes and the Index Theorem can be recovered. It turns out 
that even on a very small lattice (6x6), the exact zero modes with definite chirality are 
reproduced to a very high accuracy and the Index Theorem is satisfied exactly. For smooth 
background guage fields, the convergence of the recursion formula eq.(|4]) is indeed quite fast. 
If we require the absolute error of each matrix element satisfying the criterion 

\Y* - {DlD w ) t3 \ < e (23) 

then it usually takes less than 5 iterations to reach e = for smooth gauge configura- 
tions. For very rough configurations, the convergence could be slow and the accumulation of 
roundoff errors in successive iterations could cause the algorithm unstable. For exceptional 
configurations, det(D w ) = 0, then the inverse square root operator in D h is not defined. 
This is consistent with the fact that eq.@ breaks down when D^Dw is not positive defi- 
nite. In fact, the convergence of eq.@ has become very poor in the vicinities of exceptional 
configurations. 

After the matrix is computed, we solve the following eigenproblem 

££0f(s,i/)#(i/)=A,# , ( a ;) s = l,---,2N 1 N 2 (24) 

y /3 

where x and y are site indices, a and j3 are Dirac indices, and 4> s is normalized eigenfunction 

EE(C(«W = i (25) 

x a 

In matrix notations, these equations are rewritten as 

D h <j) s = X s cf) s (26) 
1 (27) 



Before we solve for the eigenvalues and eigenf unctions of D^, we discuss some of their general 
analytical properties in the following. Since 75^75 = D^, it follows that 

75^75 = D{ (28) 
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and the secular equation 

det(D h - Al) = det{D\ - A*!) = det[y s {D h - A*I)7 5 ] = cfet(A, - A*l) = (29) 

implies that the eigenvalues X s are either real or come in complex conjugate pairs. Using 
eq.(|28|) and eq.(|2|), we obtain 

D\D h = D h D\ D h is normal (30) 

and 

D} l + D h = D{D h (31) 

Since is normal, and D k have common eigenfunction and their eigenvalues come in 
complex conjugate pairs, 

D{<j )s = \:<p s (32) 
and the eigenvectors {4> s } form a complete orthonormal set. Eq . fl3l|) implies that 

A: + A s = A:A s ^|A s -1| 2 = 1 (33) 

Thus the eigenvalues of Dh fall on a unit circle with center at 1 and have the reflection 
symmetry with respect to the real axis. We define the chirality of an eigenmode to be 

Xs = filth = E EEtC^rTs^W (34) 

x a Q 



Using eq. (|28| ) and eq. (p2|), we obtain 

Dhl^s = A*7 5 0s (35) 
Multiplying <p\ on both sides and using the eigenvalue equation, we get 

\ S <t>\l5<t>s = K<t>ln><t>s (36) 

Then 

Xs = til^s = if A s ^ A: (37) 

The chirality of any complex eigenmode is zero. If X s is real ( or 2 ), eq.(|35|) and eq. (p6|) 
imply that 4> s has definite chirality +1 or —1. 

75& = ±& if A s = A: (38) 

It is remarkable that the zero modes and the +2 eigenmodes are both chiral. This is true for 
any D satisfying the Ginsparg- Wilson relation eq. (|2[) and the adjoint condition eq. (|28|). 
Another useful property of chirality is that total chirality of all eigenmodes must vanish , 



E^ = Y.tii^s 



s x a p 

EEE75^ = (39) 

x a p 
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where the completeness relation 

(y) = S a % y (40) 

S 

for the eigenfunctions of Dh has been used. From eq. (|39"D , we immediately obtain 

n 2 ~ n 2 = n - ~ n + (41) 

The sum of chirality of all +2 eigenmodes is equal to the index of the zero modes . Therefore 
we can identify the chirality of +2 eigenmodes as the index of the zero modes for any D 
satisfying the Ginsparg- Wilson relation and the adjoint condition. On the other hand, the 
standard Wilson-Dirac lattice fermion operator D w satisfies the adjoint condition but not 
the Ginsparg- Wilson relation. The chirality of its complex eigenmodes is zero but its real 
eigenmodes do not have definite chirality. Moreover, it does not have exact zero modes 
even in a smooth background gauge field with non-zero topologically charge. Unlike those 
D satisfying Ginsparg- Wilson relation has real eigenvalues only at or 2, D w has real 
eigenvalues at several different values and the total chirality of all ( real ) eigenmodes does 
not vanish. Therefore, strictly speaking, we should not identify the chirality of a subset of 
real eigenmodes of D w to be the index of zero modes. Recently, Gattringer, Hip and Lang 
|T8|j conjectured that the sum of the chirality of the real eigenmodes of the Wilson-Dirac 
hopping matrix in the vicinity of +2D ( D is the dimensionality ) is equal to the minus of 
the topological charge ( — Q ) of the background gauge field. However, the Wilson-Dirac 
fermion operator does not have exact zero modes without tuning the mass parameter for 
each gauge configuration. Only in the infinite volume limit, the +2D modes of Wilson-Dirac 
hopping matrix are equivalent to the zero modes of the Wilson-Dirac fermion operator with 
zero bare mass, it implies that in the infinite volume limit the approximate zero modes 
of the Wilson-Dirac operator will become exactly massless with definite chiralities and the 
Index Theorem can be satisfied exactly. 

In the following we derive some basic properties of the unitary matrix V = Dh — I 
which has eigenvalues A s — 1. The unitarity of V implies that its eigenvalues are unimodular 
( i.e., e Ws ) and the eigenvalues of Dh are in the form 1 + e Ws . The zero modes of Dh 
correspond to those of eigenvalues -1 of V. Since the eigenvalues of V are either real ( +1 
or -1 ) or come in complex conjugate pairs, it follows that 

det(V) = (_i)(™-+"+) = (_!)(«--»+) = (_i)(4-™ 2 ") (42) 

where eq. (|4T]) has been used in the second equality. Eq. (|42|) can serve as a check for the 
consistency of the eigenvalues. 

In Table 1, the eigenvalues of Dh are listed for 6 x 6 lattice in the background gauge 
field [ eqs.(|7|) and (H) ] with topological charge Q — 1; and harmonic parts with hi = 0.1 and 
hi = 0.2, and the local parts with = 0.3, = 0.4 and n\ = n 2 = 1. It is evident that 
the eigenvalues are either real ( or 2 ), or come in complex conjugate pairs in a very precise 
manner. This property is vital to obtain real and positive fermionic determinant. There is 
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one zero mode with exactly zero eigenvalue and chirality — 1. Therefore we have n_ = 1 and 
n + = and the Index Theorem Q = n_ — n + is satisfied exactly. Moreover, the Vanishing 
Theorem Q > -<=>- n + = which only holds in two dimensions is also satisfied. We also 
note that eq. ( |39"D and eq. fl4"I|) are both satisfied. All analytical properties of the eigensystem 
we have discussed above are satisfied exactly. The eigenvalues in Table 1 are plotted in Fig. 
1. 

In Table 2, the topological charge Q = —2 while other parameters are the same as in 
Table 1. There are two exact zero modes of chirality +1. Therefore we have n_ = and 
n + = 2 and the Index Theorem Q = n_ — n + and the Vanishing Theorem Q < -<=>- n_ = 
are satisfied. Again all analytical properties are satisfied exactly. The eigenvalues in Table 
2 are plotted in Fig. 2. 

We have tested many different smooth gauge configurations by changing the topological 
charge as well as other parameters in eqs. ( |7| ) and ( [8| ). The exact zero modes are always 
reproduced such that the Index Theorem and the Vanishing Theorem are satisfied. All 
analytical properties we have discussed above are satisfied exactly. We summarize some of 
our results in Table 3. 

We also investigated the robustness of zero modes and the index theorem under local 
fluctuations by varying the amplitudes Ai and A^ as well as the frequencies and 
of the background gauge field. We found that the exact zero modes are very robust under 
variations of the background. For example, for Q = 1 and n\ — n-i — 1, the exact zero modes 
and the Index Theorem are reproduced for A^ and A% ranging from 0.0 to 0.9. 

The stability of zero modes under random fluctuations are also investigated. The rough 
gauge configurations are obtained from the smooth ones by multiplying each link variable 



with a random phase [17 



U,(x) r = e ire ^U,(x) (43) 

where O^x) is a uniformly distributed random variable in (— it, n) and r controls the size 
of roughness. We found that for gauge configurations with topological charge \Q\ < 6 and 
A^ = 0.3, A^ = 0.4 and n\ = 112 = 1, the exact zero modes satisfying the Index Theorem 
can be reproduced provided that |r| < 0.3. 



4 Fermionic Determinants 

The fermionic determinant det(D) is proportional to the exponentiation of the one-loop 
effective action which is the summation of any number of external sources interacting with 
one internal fermion loop. It is one of the most crucial quantities to be examined in any 
lattice fermion formulations. The determinant of is the product of all its eigenvalues 

det(D h ) = n(l + e ies ) = (1 + e^+^dehiDh) (44) 
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Re(A) 


Im(A) 


chirality 


Re(A) 


Im(A) 


chirality 


2.00000000 


.00000000 


1.0000000 








1 99981020 


01948230 

■ V 7 1 > 7 I l. 7 KJ V 7 


00000000 


1 99981020 


- 01948230 


00000000 

a V 7 V 7 V 7 V 7 V 7 V ' V 7 V / 


1.98951740 


.14441370 


.00000000 


1.98951740 


-.14441370 


.00000000 


1.96345153 


.26788272 


.00000000 


1.96345153 


-.26788272 


.00000000 


1.93881783 


.34441410 


.00000000 


1.93881783 


-.34441410 


.00000000 


1.93303732 


.35977961 


.00000000 


1.93303732 


-.35977961 


.00000000 


1.93041388 


.36651060 


.00000000 


1.93041388 


-.36651060 


.00000000 


1.91625509 


.40059533 


.00000000 


1.91625509 


-.40059533 


.00000000 


1.90283899 


.42997879 


.00000000 


1.90283899 


-.42997879 


.00000000 


1.89600810 


.44403770 


.00000000 


1.89600810 


-.44403770 


.00000000 


1.87536396 


.48346451 


.00000000 


1.87536396 


-.48346451 


.00000000 


1.86076634 


.50900030 


.00000000 


1.86076634 


-.50900030 


.00000000 


1.84065148 


.54157648 


.00000000 


1.84065148 


-.54157648 


.00000000 


1.83027494 


.55735405 


.00000000 


1.83027494 


-.55735405 


.00000000 


1.80483165 


.59350317 


00000000 


1.80483165 


-.59350317 


.00000000 


1.78430700 


.62037289 


.00000000 


1.78430700 


-.62037289 


.00000000 


1 76077457 


64901623 

• V 7 1 ' / V 7 1 \ } 1 } 


00000000 

.V/V/V/V/V/V/V/V/ 


1 76077457 


- 64901623 


.00000000 


1.73598092 


.67700228 


.00000000 


1.73598092 


-.67700228 


.00000000 


1.71413711 


.70000585 


.00000000 


1.71413711 


-.70000585 


.00000000 


1.68715242 


.72651329 


.00000000 


1.68715242 


-.72651329 


.00000000 


1 65416649 

A . \J 7 1 1 V 7 V 7 1 ' / 


75635058 


00000000 

.V/V/V/V/V/V/V/V/ 


1 65416649 

J a V/ • 7 1 1 V / V / 1 ' / 


- 75635058 


00000000 

a V 7 V 7 V 7 V 7 V 7 V ' V 7 V 7 


1.61565168 


.78801841 


.00000000 


1.61565168 


-.78801841 


.00000000 


1.57389192 


.81893105 


.00000000 


1.57389192 


-.81893105 


.00000000 


1.51607964 


.85654060 


.00000000 


1.51607964 


-.85654060 


.00000000 


1.47365734 


.88070922 


.00000000 


1.47365734 


-.88070922 


.00000000 


1.40247961 


.91542895 


.00000000 


1.40247961 


-.91542895 


.00000000 


1.30540300 


.95222319 


.00000000 


1.30540300 


-.95222319 


.00000000 


1.26428480 


.96444468 


.00000000 


1.26428480 


-.96444468 


.00000000 


1.15140180 


.98847230 


.00000000 


1.15140180 


-.98847230 


.00000000 


.98655620 


.99990963 


.00000000 


.98655620 


-.99990963 


.00000000 


.89820891 


.99480580 


.00000000 


.89820891 


-.99480580 


.00000000 


.78417724 


.97643256 


.00000000 


.78417724 


-.97643256 


.00000000 


.65832653 


.93981873 


.00000000 


.65832653 


-.93981873 


.00000000 


.53737749 


.88655537 


.00000000 


.53737749 


-.88655537 


.00000000 


.34130700 


.75241181 


.00000000 


.34130700 


-.75241181 


.00000000 


.10715386 


.45036182 


.00000000 


.10715386 


-.45036182 


.00000000 


.00000000 


.00000000 


-1.0000000 









Table 1: The eigenvalues of Dh in a smooth background gauge field of topological charge 
Q = 1. The values of other parameters are h\ = 0.1, h 2 = 0.2, = 0.3, A 2 = 0.4 and 
rt\ — n 2 — 1. The spectrum shows that n_ = 1 , n + = and the Index Theorem and 
Vanishing Theorem are satisfied exactly. 
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Re(A) 


Im(A) 


chirality 


Rc(A) 


Im(A) 


chirality 


2.00000000 


.00000000 


-1.0000000 








2.00000000 


.00000000 


-1.0000000 








1.99989124 


.01474800 


.00000000 


1.99989124 


-.01474800 


.00000000 


1.99940122 


.03460068 


.00000000 


1.99940122 


-.03460068 


.00000000 


1.98160401 


.19092820 


.00000000 


1.98160401 


-.19092820 


.00000000 


1.94402467 


.32987486 


.00000000 


1.94402467 


-.32987486 


.00000000 


1.93390001 


.35753429 


.00000000 


1.93390001 


-.35753429 


.00000000 


1.91210916 


.40994742 


.00000000 


1.91210916 


-.40994742 


.00000000 


1.90373008 


.42810273 


.00000000 


1.90373008 


-.42810273 


.00000000 


1.90165787 


.43245009 


.00000000 


1.90165787 


-.43245009 


.00000000 


1.89104063 


.45392355 


.00000000 


1.89104063 


-.45392355 


.00000000 


1.87807847 


.47851667 


.00000000 


1.87807847 


-.47851667 


.00000000 


1.85325494 


.52149401 


.00000000 


1.85325494 


-.52149401 


.00000000 


1.83821525 


.54533952 


.00000000 


1.83821525 


-.54533952 


.00000000 


1.81009961 


.58629226 


.00000000 


1.81009961 


-.58629226 


.00000000 


1.79786888 


.60283103 


.00000000 


1.79786888 


-.60283103 


.00000000 


1.76948744 


.63866194 


.00000000 


1.76948744 


-.63866194 


.00000000 


1.74331258 


.66894425 


.00000000 


1.74331258 


-.66894425 


.00000000 


1.71748554 


.69657340 


.00000000 


1.71748554 


-.69657340 


.00000000 


1.69419796 


.71978413 


.00000000 


1.69419796 


-.71978413 


.00000000 


1.66871002 


.74352331 


.00000000 


1.66871002 


-.74352331 


.00000000 


1.63944000 


.76884100 


.00000000 


1.63944000 


-.76884100 


.00000000 


1.59262521 


.80547834 


.00000000 


1.59262521 


-.80547834 


.00000000 


1.54820232 


.83634575 


.00000000 


1.54820232 


-.83634575 


.00000000 


1.48363430 


.87527017 


.00000000 


1.48363430 


-.87527017 


.00000000 


1.43095251 


.90237461 


.00000000 


1.43095251 


-.90237461 


.00000000 


1.39959189 


.91669315 


.00000000 


1.39959189 


-.91669315 


.00000000 


1.26607166 


.96395325 


.00000000 


1.26607166 


-.96395325 


.00000000 


1.18891548 


.98199335 


.00000000 


1.18891548 


-.98199335 


.00000000 


1.06804551 


.99768222 


.00000000 


1.06804551 


-.99768222 


.00000000 


.90949855 


.99589632 


.00000000 


.90949855 


-.99589632 


.00000000 


.87904797 


.99265835 


.00000000 


.87904797 


-.99265835 


.00000000 


.72327655 


.96094960 


.00000000 


.72327655 


-.96094960 


.00000000 


.56601895 


.90092200 


.00000000 


.56601895 


-.90092200 


.00000000 


.48688458 


.85831962 


.00000000 


.48688458 


-.85831962 


.00000000 


.18512327 


.57963429 


.00000000 


.18512327 


-.57963429 


.00000000 


.00000000 


.00000000 


1.0000000 








.00000000 


.00000000 


1.0000000 









Table 2: The eigenvalues of Dh in a background gauge field of topological charge Q = -2. 
The values of other parameters of the gauge field are the same as in Table 1. The spectrum 
shows that n_ = 0, n + = 2 and the Index Theorem and Vanishing Theorem are satisfied 
exactly. 
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Q 


n + 


n_ 


_L 

rq 


712 


-5 


5 








5 


-4 


4 








4 


-3 


3 








3 


-2 


2 








2 


-1 


1 








1 

















1 





1 


1 





2 





2 


2 





3 





3 


3 





4 





4 


4 





5 





5 


5 






Table 3: The zero modes versus the topological charge. The Index Theorem and the Van- 
ishing Theorem are satisfied exactly. Eq.(^TJ) is also satisfied. 

where deti(D h ) is equal to the product of all non-zero eigenvalues. Since the eigenvalues 
are either real of come in complex conjugate pairs, det(D h ) must be real and positive. For 
Q = 0, then n + + ri_ = and det(D h ) = deti(D h ). For Q ^ 0, then n + + n_ ^ and 
det(Dh) = 0, but deti(Dh) still provides important information about the spectrum. In 
continuum, exact solutions of fermionic determinants in the general background U(l) gauge 
fields on a torus was obtained by Sachs and Wipf |H|. In the following we compute deti(Dh) 
for several different gauge configurations and then compare them with the exact solutions 
in continuum. For simplicity, we turn off the harmonic part ( hi = h 2 = 0.0 ) and the local 
sinusoidal fluctuations ( = = ) in Eqs. (0) and ([|) and examine the change of 
deti(Dh) with respect to the topological charge Q. For such gauge configurations, the exact 
solution [[H]] is 

det 1 [D(Q)]=N ] J(^) lQl (45) 
where the normalization constant iV is fixed by 




such that deti[D{\)] = 1. 

In Table 4, the fermionic determinants deti(Dh) are listed for 8x8 and 16 x 16 lattices 
respectively. They agree with the continuum exact solutions very well for small Q but the 
error goes up to 10% for large Q ~ 10. For a fixed Q, the error decreases with respect to 
the increasing of the size of the lattice. 

We have also computed the fermionic determinant deti[Dh] for more general gauge 
configurations with non-zero topological charge, harmonic parts and local parts; as well as 
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8x8 


16 x 16 


Q 


deti[D(Q)]exact 


deti[D h (Q)] 


deti[D(Q)] exact 


deti[D h (Q)) 


1 


1.00000 


1.00000 


1.00000 


1.00000 


2 


2.82843 


2.77348 


5.65685 


5.66186 


3 


6.15840 


5.96891 


24.6336 


24.0615 


4 


11.3137 


10.7157 


90.5097 


90.4894 


5 


18.3179 


16.9340 


293.086 


286.003 


6 


26.8177 


24.2001 


858.166 


822.664 


7 


36.1083 


32.0006 


2310.93 


2170.94 


8 


45.2548 


40.2920 


5792.62 


5354.28 


9 


53.2732 


45.7353 


13637.9 


12309.5 


10 


59.3164 


50.2816 


30370.0 


27336.6 



Table 4: The fermionic determinant versus the topological charge Q. The normalization 
constant is chosen such that deti[Dh(l)} = 1. The results for the 8x8 lattice are listed in 
the second and the third columns, while for the 16 x 16 lattice in the last two columns. The 
exact solutions are computed according to Eq.( ^ ). 

in the topologically trivial sector with only harmonic parts or local parts. They are all in 
very good agreement with the continuum exact solutions [[Uj. 



5 Conclusions and Discussions 

In this paper, we have implemented the square root operator in Neuberger's proposal [l[] of 
exactly massless quarks on the lattice by the recursion formula eq.(^). Our numerical tests 
in two dimensional background gauge fields assert that D h indeed reproduces the exactly 
massless fermions on a torus. In particular, for smooth background gauge fields with non- 
zero topological charge, the exact zero modes with definite chirality are reproduced and the 
Index Theorem is satisfied exactly. 



Due to the recent work of Neuberger |l|, |2|, Hasenfratz, Laliena and Niedermayer |TT 



and Liischer ||12|| , we now have a better understanding of the chiral symmetry on the lattice. 



The Nielson-Ninomiya theorem [13| can be circumvented by replacing the continuum chiral 
symmetry {D, 75} = by the Ginsparg- Wilson relation {D, 75} = -D75D on the lattice. The 
Ginsparg- Wilson relation gaurantees that the correct continuum anomaly can be recovered 
on the lattice. Currently, there are two classes of solutions satisfying the Ginsparg- Wilson 
relation. One of them is the explicit solution obtained by Neuberger [Q, which grew 
out of the overlap formalism || |], ^|, |J. If one attempts to solve the Ginsparg- Wilson 
relation by writing D = 1 + T, one must end up with the solution (75T) 2 = 1. Naively 
one can choose any lattice operator T which satisfies this condition and gives the correct 
continuum action in the continuum limit. The choice of T = D^D^D^ 1 / 2 where D w is the 
standard Wilson-Dirac action with negative mass term is identical to Neuberger's overlap 
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solution. We note in passing that for any non-trivial T, the inverse square root operation 
in T seems to be inevitable . Another class of solutions of Ginsparg- Wilson relation is 
the fixed point lattice Dirac operator h FP ]19 which is the implicit solution of the non- 



linear classical saddle-point equations, and is obtained by recursive iterations. The technical 
difficulty of that approach is how to parametrize h FP such that it is sufficiently precise but 
the computational costs are still affordable. It is instructive to compare our numerical results 
on zero modes with those obtained by Farchioni and Laliena |2(J] using fixed point lattice 
Dirac operator. It is evident that their zero modes are not exactly zero. The discrepancies 
are essentially due to the truncation errors in the parametrization of h FP and the round- 
off errors in the computations performed at each step of the iterations. These errors are 
intrinsically mixed together and therefore are very difficult to control. Furthermore, we 
suspect that the fermionic determinants deti(h FP ) of their h FP would have relatively large 
errors. On the other hand, we do not have any serious technical difficulties to compute Dh 
to a very high precision using the recursion formula eq. (f|). Therefore we can obtain exact 
zero modes and fairly accurate fermionic determinants even on a finite lattice. 

It is interesting to note that tracing the Index Theorem on the lattice has a long history. 



In 1987, Smit and Vink [[T7| investigated the spectrum of Wilson and Staggered fermions in 
topologically non-trivial background gauge fields in two dimensions. Some remnants of the 
index theorem were found but exact zero modes with definite chirality were not obtained. 
The Index Theorem on the lattice was also recognized by Neuberger and Narayanan || [H| 
in the context of overlap formalism. The index was computed by studying the level crossings 
in the spectral flow of H(m) = r y^D w {m) as a function of m and was verified numerically to 
equal to the topological charge of the background gauge field. When the lattice Dirac fermion 
operator D h was proposed by Neuberger |], it became clear that the index computed from 
level crossings of H(m) is equal to the index of Dh- 

Although the recursion formula eq.(f|) is quite effective in computing Dh for exactly 
massless fermions in two dimensions, however, we expect that it would still be expensive 
for four dimensional lattice QCD simulations. An algorithm which invloves at most matrix 
multiplications but without matrix inversions at every step of iterations would be the most 



desirable. We are now contemplating such possibilities [22]. There are three emerging 
investigations that we would like to carry out in the near future. The first is the spectrum 
of Dh in the four dimensional topologically non-trivial background gauge fields, which is 



essentially an extension of the present investigation to four dimensions j23]. The second is 
the quenched QCD calculations of some chirally sensitive observables such as the kaon weak 
matrix elements which were also calculated by Blum and Soni |24fl using the Domain Wall 



quarks |25| , j26| . Since the Domain Wall quarks in the limit of an infinite fifth dimension 
becomes exactly massless and is equivalent to D h , it would be interesting to compare future 



results of using Dh with those given by Domain Wall quarks with finite fifth dimension p4 
and examine whether any improvements can be made by using Dh with exactly massless 
quarks. The third is to perform dynamical fermion simulations of Dh in two dimensions, for 
example, the Schwinger model. If the results continue to agree with the continuum exact 
solutions, we would attempt to confront one of the most challenging problems in QCD such 
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as the chiral symmetry breaking by measuring the spectral density of the eigenmodes of 
near zero but not exactly zero. 
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Figure 1: The eigenvalues of D h in Table 1. 
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Figure 2: The eigenvalues of D h in Table 2. 
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